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The results on chaos in FRW cosmology with a massive scalar field are extended to 
another scalar field potential. It is shown that for sufficiently steep potentials the 
chaos disappears. A simple and rather accurate analytical criterion for the chaos 
to disappear is given. On the contrary, for gently sloping potentials the transition 
to a strong chaotic regime can occur. Two examples, concerning asymptotically 
flat and Damour-Mukhanov potentials are given. 

1 Introduction 

The studies of chaotical dynamics ofxlosed isotropic cosmological model has a long 
story. They were initiated by papers u where the possibility to avoid a singularity at 
the contraction stage in such a model with a minimally coupled massive scalar field 
was discovered. Later it was found that this model allows the existence of periodical 
trajectories □ and aperiodical infinitely bouncing trajectories having a fractal nature 
S. In □ the set of periodical trajectories was studied from the viewpoint of dynamical 
chaos theory. It was proved that the dynamics of a closed universe with a massive 
scalar field is chaotic and an important invariant of the chaos, the topological 
entropy, was calculated. 

However, most of modern scenarios based on ideas of the string theory and 
compactification naturally lead to another forms of potential which are exponential 
or behave as exponential for large ip (see, for example, the paper of Gunther and 
Zhuk □). This steepness of the potential apparently changes the possibilities of 
escaping the singularities and alters the structure of infinitely bouncing trajectories. 
Under some conditions, which are described in Section 2, the chaotic behaviour can 
completely disappear. 

One the other hand, potential, less steep than quadratic-, one can give rise a 
chaotic dynamics which differ qualitatively from described in □. Two examples are 
presented in Section 3. 

2 Chaotic properties of closed FRW model with a scalar field 

We shall consider a cosmological model with an action 

S = Jd*xy/^{&R+±<rd„<pd v <p-V(<p)y (1) 

For a closed Friedmann model with the metric 

ds 2 =dt 2 -a 2 (t)d 2 n {3 \ (2) 
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where a(t) is a cosmological radius, cPtl^ is the metric of a unit 3-sphere and with 
homogeneous scalar field ip we have the following ODE system: 

™k(a+*+±)+!!£ -^M = 0, (3) 
16i\ 2a 2oj 8 4 ' w 

+ ^ +V '( ip ) = O , ( 4 ) 
a 

where the scalar field potential V((f)) is a smooth nonnegative function with V(0) = 
0. 

This system has one first integral of motion 



3 



-m 



+ + + 2V{ V )) = o. (5) 

It is easy to see from Eq. (5) that the points of maximal expansion and those 
of minimal contraction, i.e. the points, where a = can exist only in the region 
where 

Sometimes, the region defined by inequalities (6) is called Euclidean, and the op- 
posite region is called Lorentzian. We will useJbhis definition for brevity, though 
we consider these notations rather meaningless U. 

It also can be shown that the possible points of maximal expansion are localized 
inside the region 

while the possible points of minimal contraction lie outside this region (7) being at 
the same time inside the Euclidean region (6) 0. 

In our paper u it was shown that the region of possible points of maximal 
expansion has quite a regular structure. In the left, closest to axis a — part of 
this region there are points of maximal expansion after which trajectory goes to 
singularity. Then one can see the region where after the going through the point 
of maximal expansion a trajectory undergoes the "bounce" i.e. goes through the 
point of minimal contraction. Then we have the region where after going through 
the point of maximal expansion trajectory has a "tp - turn " i.e. has the extremum 
in the value of the scalar field tp and then falls into a singularity. Then one has 
the region corresponding to the trajectories having bounce after one oscillation in 
(p and so on. 

To avoid a misunderstanding, let us indicate once more what initial condition 
space we use. When we start from maximal expansion point, we fix one time 
derivative (d = 0). So, for fix the initial condition completely, we need only to 
specify initial values of a and (p. The initial (p is determined from the constraint 
equation (5), and our initial condition space becomes (a, (p). We will study the 
structure of this space, investigating the location of points of maximal expansion, 
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starting from which a trajectory has a bounce, but not the location of bounce itself. 
The ip = cross-section of regions, leading to bounce, are called as "intervals". 

Remembering that a trajectory describing an expanding universe must have a 
point of maximal expansion, we can further apply this analysis to bouncing tra- 
jectories. Their second point of maximal expansion may lie either inside "bounce" 
regions, defined at the first step or between them. This fact generates a more tiny 
substructure of the region under consideration. This substrucure of regions having 
two bounces repeats the general structure of regions having at least one bounce 
and so on and so forth. Continuing this process ad infinitum we can get the fractal 
zero-measure set of infinitely bouncing trajectories escaping the singularity. 

Numerical investigations show also that all simple periodical trajectories (i.e. 
having only one bounce per period) have a full stop point on the Euclidean boundary 
(the curve defined by the equality in (6)). Moreover, trajectories, going from the 
boundary inside the Euclidean region has a point of maximal expansion almost 
immediately and then go towards a singularity. So, periodical trajectory approaches 
their bounce point on the boundary from the Lorentzian side. Hence, we have a 
necessary condition for given region of the Euclidean boundary to contain a full 
stop points of periodical trajectories (it is the condition that a trajectory starting 
with zero velocities from the Euclidean boundary goes into the Lorentzian region) : 

% < (8) 
a da 

where function tp(a) in the right-hand side is the equation of the Euclidean boundary 
(6). 

The point on the Euclidean boundary where a trajectory starting with zero 
velocities have a direction tangent to the Euclidean boundary can be obtained 
substituting the equality in (6) into the equations of motion. It was first introduced 
by Page El for massive scalar field potential V{4>) = m 2 (j) 2 /2. In this case 




Ppage — 

dpage = 1/™, (9) 

except for the trivial solution ip = 0, a — oo. 

The first periodical trajectory have full stop point at ip = 1.19(p page and 
a = 0.83a page . These obtained numerically values, rather close to the Page point, 
represent a right-down boundary of a set of full stop points of periodical trajectories 
on the Euclidean boundary. 

Using the equation of motion (3)-(4) the criterion (8) gives for an arbitrary 
scalar field potential 



V(<p)>yj^V'(<p) (10) 

The condition (10) may be treated as restricting a local steepness of the function 
V (</)). In can be easily seen that for power-law function the Page value 4> page such 
that for all <f) > <j) pag e the condition (8) is satisfied, always exists. For steeper poten- 
tials the situation changes. For example, the potential V(ip) = Mq (cosh(ip / (po) — l)i 
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Figure 1. Example of trajectories with the initial conditions close to the boundary separating 
trajectories falling into ip = +00 (trajectory 1) and ip = —00 (trajectories 2 — 5) singularities 
for the case <po < ^J- mp. This boundary is sharp, no fractal structure is present. Trajectories 
2 — 5 have a zigzag-like form, no periodical trajectories are present. The long-dashed line is the 
Euclidean boundary, the short-dashed line is the separating curve. 



studied in □ has a Page point only if ip > -^=mp. In the opposite case the con- 
dition (8) is never satisfied and all the trajectories starting from the Euclidean 
boundary go into the Euclidean region and soon experience the point of maximal 
expansion. It was confirmed numerically that the chaos is absent in this case (see 
Fig.l.). 

For steeper potentials like 

V{fp) = M 4 (exp(^ 2 /Vo) + cxp( _ ¥ ,2 / ^2 ) _ 2) 

inequality (8) is definitely violated for large <j) but, depending on <j)Q it can be 
satisfied for intermediate <j). For different 0o there are two or zero Page points in this 
case. The value of 4>q corresponding to the Page points disappearance (0.905mp) 
differ from the 4>o of the chaos disappearence obtained numerically (0.96mp) for 
~ 6% 13. It is interesting that when the chaos exists, the number of bounce intervals 
is finite and the whole picture of the chaos is similar to the picture obtained for a 
system with a massive scalar field and a hydrodynamical matter E£l. 

So, studying the possibility to satisfy the condition (8) at the Euclidean bound- 
ary, we have an easily calculated and rather accurate criterion for the existence of 
the chaotic dynamics in the system (3)-(5). 
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3 Chaotic dynamics for gently sloping potentials 



In this section we describe the opposite case - the potential which is less steep than 
the quadratic one. We will see that in this case the transition to a qualitatively 
stronger chaos may occur. 

3.1 Asymptotically flat potentials and merging of the bounce intervals 
Let us start with the potential 

2 

%) = M 4 (l-exp(^)), (11) 

where M and ipo are parameters. M determines the asymptotical value of the 
potential for ip — ► ±oo. 

It can be easily checked from the equations of motion that multiplying the 
potential to a constant (i.e. changing the Mo) leads only to rescaling a. So, this 
procedure do not change the chaotic properties of our dynamical system. On the 
contrary, this system appear to be very sensitive to the value of <po. We plotted in 
Fig. 2. the <p = cross-section of bounce intervals depending on <p a . In our numerical 
investigations we use the units in which mp/y/T&K = 1, and below we present our 
results in these units because now most of the interesting events occur for the range 
of parameters of the order of unity. This plot represents a situation, qualitatively 
different from studied previously for potentials like V ~ ip 2 and steeper. Namely, 
the bounce intervals can merge. 

Let us see more precisely what does it means. For <po > 0.82 the picture is 
qualitatively as for a massive scalar field - trajectories from 1-st interval have a 
bounce with no <p-turns before it, trajectories which have initial point of maximal 
expansion between 1-st and 2-nd intervals fall into a singularity after one <^-turn, 
those from 2-nd interval have a bounce after one y-turn and so on. For ipo a bit 
smaller than the first merging value the 2-nd interval contains trajectories with two 
y-turns before bounce, the space between 1-st interval (which is now the product of 
two merged intervals) and the 2-nd one contains trajectories falling into a singularity 
after two c^-turns. There arc no trajectories going to a singularity with exactly one 
yj-turn. Trajectories from the 1-st interval can experience now a complicated chaotic 
behavior which can not be described in as similar way as above. 

With (po decreasing further, the process of interval merging being to continue 
leading to growing chaotisation of trajectories. When n intervals merged together, 
only trajectories with at least n oscillations of the scalar field before falling into 
a singularity are possible. Those having exactly n 92-turns have their initial point 
of maximal expansion between 1-st bounce interval and the 2-nd one (the second 
interval now contains trajectories having a bounce after n 93-turns). For initial 
values of the scale factor larger then those from the 2-nd interval, the regular 
quasiperiodic structure described above is restored. 

Numerical analysis shows also that the fraction of very chaotic trajectories as a 
function of ipo grows rapidly with tpo decreasing below the first merging value. To 
illustrate this point we studied the behavior of trajectories starting from the point 
of maximal expansion with initial values of a located in the range of the first two 
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Figure 2. The (p = cross-section of the bounce intervals for the potential (11) depending on 
ipo- Consecutive merging of 5 first intervals can be seen in this range of ipg. 

intervals on the greed with step 0.002. Total number of trajectories for each ipo is 
equal to 1000. We plotted in Fig. 3 the number of trajectories which do not fall into 
a singularity during first 50 oscillations of the scalar field <p. We do not include 
trajectories with the next point of maximal expansion located outside the 2-nd 
(or the 1-st one, if merging occurred) interval, so all counted trajectories avoid a 
singularity during this sufficiently long time interval due to their extreme chaoticity, 
but not due to reaching the slow-roll regime. Before merging, the measure of so 
chaotic trajectories is extremely low and they are undistinguishablc on our grid. 
When ip becomes slightly lower than the value of the first merging, this number 
begin to grow rather rapidly and for (po ~ 0.6 near 10% of trajectories from the 1-st 
interval on our grid experience at least 50 oscillation before falling into a singularity. 

We point out that for a simple massive scalar field potential only <~ 10~ 2 trajec- 
tories in the same range of the initial scale factors have at least one bounce (using 
the abovementioned quasiperiodicity in the structure of initial condition space (a, 
ip) we can estimate this measure as a ratio of interval width and distance between 
near-by intervals). The width of subintervals containing trajectories not falling into 
a singularity after only one bounce is about one hundred times less than the width 
of " main" intervals and so on. The common numerical calculation accuracy is un- 
sufficient for distinguishing even the sole trajectory with 50 oscillations and a being 
in the range of first two intervals. 

In contrast to this, the chaos for the potential (11) is really significant. Detail 
of intervals merging including the description out of <p> = cross-section require 
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Figure 3. Number N of trajectories do not falling into a singularity during 50 oscillating times 
for the potential (11) depending on the parameter tpo. The scale factor of the initial maximal 
expansion point varies in the range of the 1-st and 2-nd intervals which merge at ipo = 0.82. Total 
number of trajectories is equal to 1000. 



further analysis. 

For large initial a the configuration of bounce intervals for potential (11) looks 
like the configuration for a massive scalar field potential with the effective mass 
rrieff = (V2AIq)/(Pq. The periods of corresponding structures coincides with a 
good accuracy (mAa ~ 2.8 in the initial condition space) though the widths of the 
intervals for the potential (11) is bigger then for V = (m 2 e ^ip 2 )/2. 

3.2 D amour- Mukhanov potentials 

The very chaotic regime described above is possible also for potentials, which are not 
asymptotically flat, if the potential growth is slow enough. We will illustrate this 
point describing a particular (but rather wide) family of potentials having power- 
low behavior - Damour - Mukhanov potentials □. They was originally introduced 
to show a possibility to have an inflation_behaviour without slow-roll regime After, 
various issues on inflationary dynamics Ej and growth of perturbation ExEj for this 
kind of scalar field potential was studied. 

The explicit form of Damour-Mukhanov potential is 



with three parameters -Mo, q and ipo. 

For <p <C ifio the potential looks like the massive one with the effective mass 
m e ff = Mq/lpq. In the opposite case of large ip it grows like ip q . 

As in the previous section, the chaotic behavior does not depend on Mq. So, 
we have a two-parameter family of potentials with different chaotic properties. Nu- 
merical studies with respect to possibility of bounce intervals merging shows the 
following picture (see Fig. 4): for a rather wide range of q there exists a correspond- 
ing critical value of ipo such that for <po less than critical, the very chaotic regime 
exists. Increasing q corresponds to decreasing critical (fQ. 
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Figure 4. The value ^jq of the potential (12) corresponding to the first merging of the bounce 
intervals depending on q. 



Surely, since this regime is absent for quadratic and more steep potentials, q 
must at least be less than 2. We can see clearly the very chaotic regime for q < 1.24. 
The case q = 1.24 lead to strong chaos for tpo < 1-4 x 10~ 5 and the critical (fo 
decreases with increasing q very sharply at this point. We did not investigated 
further these extremely small values of ipo, because the physical meaning of such 
kind of potential is very doubtful. 
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